Inflatonic solitons in running mass inflation 



Kari Enqvist a ' b , Shinta Kasuya b , and Anupam Mazumdar c 
a Department of Physical Sciences, P. 0. Box 64, FIN-00014, University of Helsinki, Finland. 
b Helsinki Institute of Physics, P. O. Box 64, FIN-00014, University of Helsinki, Finland. 
c The Abdus Salam International Centre for Theoretical Physics, 1-34100, Trieste, Italy. 

(June 26, 2002) 

The inflaton condensate associated with a global symmetry can fragment into quasistable Q balls, 
provided the inflaton oscillations give rise to an effective equation of state with negative pressure. 
We study chaotic inflation with a running inflaton mass and show that, depending on the sign of 
the radiative mass correction, the process of fragmentation into inflatonic Q balls can actually take 
place even though no net charge exists. If the main decay channel of the Q ball is to fermions, the 
universe will be reheated slowly via surface evaporation. 
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I. INTRODUCTION 



\Q ' Reheating of the universe, in which the inflaton decay product provide a thermal bath of radiation, is one of the 
cornerstones of any inflationary model. At the end of inflation the coherent inflaton condensate may decay either 
perturbatively Q|, one quantum at the time, or non-perturbatively Q via a parametric resonance in a collective process 
denoted as preheating. Usually in the literature reheating has been considered to be a volume effect so that it tends 
to be relatively fast, and as a consequence, the reheat temperature is relatively high. In supersymmetric models this 
^ I has turned out to be somewhat problematic because a large reheat temperature generates a large gravitino abundance 

■ which is in conflict with the successful big bang nucleosynthesis j^] . 

However, recently we have pointed out [Q a novel possibility in which reheating could be realized as a surface effect, 
noting that the evaporation rate from a surface is always suppressed by the area of the evaporating surface. Therefore 
the larger the surface, the smaller the reheat temperature of the universe. Such a situation can be achieved if the 
zero-mode inflaton condensate does not remain spatially coherent but breaks down into quasi-stable lumps. This can 
r"| • happen if the inflaton field carries a global U(l) charge which protects the lumps; for instance, the inflaton could be a 
Oh! complex rather than a real scalar field. The ground state of such a condensate lump would be a non-topological soliton 
called Q ball ||, the properties of which depend on the exact form of the scalar potential. A general requirement for 
I its existence is that the scalar potential should grow slower than |$| 2 . 

Inflaton condensate fragmentation has recently been discussed in the context of hybrid inflation models both 
numerically |6| and analytically [Q. However, in these cases the shape of the potential does not render the condensate 
• i-h , lumps quasi-stable and their appearance is merely a part of the complicated preheating dynamics. In fact, there they 

■ appear to make reheating even more effective rather than slowing it down Q| . 

Q ball formation has been considered in detail in the context of the MSSM flat direction condensates (Affleck-Dine 
condensates ||), and can be either completely stable fHl| if SUSY is broken by gauge mediation, or quasi-stable 



-|16| if SUSY is broken by gravity mediation. In the Affleck-Dine case the condensate is charged up by a helical 
motion in the field space, and fragmentation takes place when the fluctuations in the real and the imaginary directions 
grow to become of order one. 

In this paper we consider the simple case of the chaotic inflation but with a complex inflaton field (which is 
rather more natural) and a running inflaton mass, such as might be generated radiatively by the Yukawa and/or 
gauge couplings of the inflaton to other fields. We demonstrate numerically that at the end of inflation the inflaton 
condensate fragments into lumps and forms Q balls. The initial fragmentation occurs because of the negative pressure 
the condensate feels during the inflaton oscillations. We show that the fluctuations in the imaginary direction grow 
because of the mode-mode coupling between the real and the imaginary parts and thus provide the motion in field 
space necessary for solitogenesis. 

The paper is organized as follows. In Section II we describe Q balls in a running mass potential and discuss how the 
fluctuations grow for a certain running mass inflaton potential in Section III. In Section IV we study Q-ball formation 
due to mode-mode couplings by lattice simulations. Section V is devoted to the implications of Q-ball formation and 
a discussion on the reheat temperature. In Section VI we conclude with a brief discussion highlighting the results and 
describing future prospects. 
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II. Q BALLS IN THE RUNNING MASS INFLATON POTENTIAL 



If the inflaton has couplings to other fields, then, in general, the inflaton mass should receive radiative corrections 
7J . The simplest chaotic inflation potential with a running inflaton mass could then be written as 



V = m 2 |$| 2 



1 + K log 



where the coefficient K could be negative or positive, and m is the bare mass of the inflaton. The logarithmic 
correction to the mass of the inflaton is something one would expect to arise because of the possible Yukawa and/or 
gauge couplings to other fields. 

Since we assume that \K\ <C 1, the second term in the bracket does not affect the inflaton dynamics or generation 
of perturbations much so that with a running mass the universe would experience the same chaotic inflationary stage 
as with the usual m 2 |$| 2 potential. However, as we shall see, the values K < lead to interesting cosmological 
consequences. Though it is not pertinent, we note that the potential Eq. (jl]) can be generated in a supersymmetric 
theory if the inflaton has a gauge coupling jl5],Q. This is due to the fact that the gaugino loops contribute as a 
negative correction which dominates. In this case the value of K is given by 

K—f^, (2) 

where mi/2 is the gaugino mass and m,~ t denotes the mass of slepton. It is also possible to obtain the potential Eq. (jl]) 
in a non-supersymmetric (or in broken SUSY) theory, provided the fermions live in a larger representation than the 
bosons. In this latter situation the value of K is determined by the Yukawa coupling pj 

h 2 

K = -C—. (3) 

Here C is the total number of fermionic loops, and h is the Yukawa coupling. As long as \K\ <C 1, the dominant 
contribution to the potential comes from m 2 |<I>| 2 term during inflation, and inflationary slow roll conditions are 
satisfied as in the case of the standard chaotic model. The generation of the amplitude of the density perturbations 
are constrained by the COBE normalization, which implies 

5p m K 
p Mp 

which results in m ~ 10 13 GeV. 

The most important virtue of the inflationary potential Eq. (|l|) is that it respects a global £7(1) symmetry, and for 
a negative K, the potential is shallower than m 2 !^ 2 . It has been well known that such a potential admits a Q-ball 
solution |ll| . This means that the energy of a Q-ball configuration is less than a collection of free scalars carrying an 
equivalent charge. The field configuration is given by t) = (e luJt (j)(x))/V2 while its energy and charge read 



E= I d 3 ^ 

Q = uj [ tfdPx . (5) 



Here we assume that the potential V(<fi) has a global minimum at the origin and is invariant under a global £7(1)- 
transformation. 

Depending on the slope of the potential it is possible to have an energetically favorable state in the form of Q balls. 
The energy of the Q ball is E = [iQ, where /i lies in the range /to < ^ < mo, and [Iq and mo are given by 



(1q = mm 



, m = y/V"(0). (6) 



The Q ball is stable against decay into other particles if fiq < rridecay, where nidecay is the smallest mass of the 
particles which shares the same £7(1) charge as the complex field and q is the charge carried by the particle. In 
this case, a Q ball exists forever. On the other hand, if the stability condition does not hold, Q balls will eventually 
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decay into lighter particles. However, if the (final) decay products are pairs of fermions, the decay should proceed 
by the evaporation of charges through the surface of Q balls because of the Pauli blocking JlJ]. Thus, inflatonic Q 
balls with a large enough charge could live very long, in which case reheating may take place very slowly, as we shall 
discuss. 

The potential of the form Eq. ([!]) generically forms a thick-wall Q balls with a Gaussian profile p| 

4>(r) oc exp(— \K\m 2 r 2 ) , (7) 
so that the size of such a Q ball can be estimated by 

R~ IKl'^m- 1 . (8) 
The energy of a Q ball can be approximately given by 

E ~ mQ , (9) 
and the total charge accumulated by a Q ball is given by 

Air 

Q ~ ui(f>lV ~ —mfilR 3 , (10) 
where (f>o is the field value at which the Q ball forms, and V ~ 4irR 3 /3 is the volume of the Q ball. 

III. AMPLIFICATION OF FLUCTUATIONS 

In the context of Affleck-Dine baryogenesis, the homogeneous mode is rotating around the origin in the potential. 
In such a case fluctuations tend to grow when on the average there is a negative pressure, which can be achieved in 
the potential Eq. (|l]) when K < |l9|,^0| . Expanding the potential in Eq. (Q) one finds 

v^^i^^y cx^*. at) 

where we assume \K\ <§; 1, since this term is actually a correction to a tree level potential. The equation of state for 
a field rotating in such a potential is then given by [] 

P- — ttP' ( 12 ) 



2 + K 



where p and p is a pressure and energy density of the scalar field, respectively. Therefore, negative value of K is 
crucial for negative pressure, which is the core of our discussion. Notice the difference between this case and the usual 
chaotic potential V = to 2 |$| 2 , where on the average the oscillating inflaton condensate acts as a pressureless fluid. 

When the homogeneous mode is rotating in the potential, it is appropriate to write a complex field as <& = {4>e l9 )/ y/2, 
and decompose it into a homogeneous part and fluctuations: <f> — > <f> + 8<j> and 9^9 + 69. The equations of motion 
of the field read as 

4> + 3Hj)-9 2 (j) + V'((l)) =0, (13) 
4>9 + 3H<f>9 + 2^)9 = 0, (14) 



for the homogeneous mode, and 



. . . V 2 

3H6(/) - 29059 - 9 2 5<p T 6(f> + V"(<j>)8<f> = 0, (15) 

a z 

3H<f>59 + 2{<j)69 + 954>) - 2^96^ - (/>\s9 = 0, (16) 



1 Because in the Affleck-Dine case the exact motion of the field is helical and hence strictly speaking non-periodic, the equation 
of state actually oscillates . 
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for the fluctuations, and 



v'{4>) = m 2 

V"{4>) = m 2 



l + K + K\og 



1 + 3K + K log 



( <t> 2 



\2M 2 
\2M 2 J 



(17) 
(18) 



Equation (|TJ) represents the conservation of the charge (or number) within the physical volume: 9<fi 2 a 3 = const. Since 
the energy density of the scalar field dominates the universe, the homogeneous part of the field evolves as 



a(i) 
a 



-3/(2+K) 



-6K/(2+K) 



m 2 . 



We are now going to find the most amplified mode. To this end, we take the solutions in the form 

' a(t) N-3/(3+*) 



ao 

a(f) 
a 



D e 



a{t)+ikx 



-3K/(2+K) 



89 e a{t)+tkx . 



(19) 
(20) 

(21) 
(22) 



If a is real and positive, these fluctuations grow exponentially, and go nonlinear to form Q balls. Putting these forms 
into Eqs.(|l5|) and (|l6|), we get the following condition for non-trivial 5<f)o and S9q, 



:iK : Ha + a + a 2 + K + 2m 2 K a - 6K ^ 2+K ^ 



2 + K 



2ma 

4>o 



~2ma- 6K /( 2+K Uo 
, k 2 3K 



3K 
2 + K' 



■H + a 



a + a 



K 



(4 - ZK)H 2 Ha 



= 0, 



(23) 



where we set ao = 1. 

Let us neglect the cosmological expansion and assume a ^ a 2 for simplicity. Then, Eq.(E3J) will be simplified as 



— + 2m 2 K -2m(poa 
a 2 

2ma n k 2 



(24) 



In order for a to be real and positive, we must have 



-=■ + 2m 2 K < 0. 



(25) 



As we are considering negative K , an instability band exists. This is because the rotating inflaton field in the potential 
flatter than 4> 2 has a negative pressure. Thus the instability band should be found in the range 



k 2 k 2 
< — ^ < = 2m 2 \K\ 



(26) 



We can easily derive that the most amplified mode lies at about the center of the band: (k most /a) 2 ~ m 2 \K\(l — \K\/A) : 
and the maximum growth rate is a mos t = |A"|m 3 / 2 /2. Q 



2 The coefficients of these results are slightly different from the previous estimations Jl^,|l3| because these were based on the 
approximation 2 ~ m 2 . Replacing this with the more accurate 8 2 = V'(4>)/4> = m 2 [l + K + K log(<j!> 2 /2M 2 )], and still assuming 
a matter-dominated universe, one finds the same estimates as in the present paper. 
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After inflation the inflaton field is just oscillating around the origin of its potential. It is thus more appropriate 
to decompose the inflaton field into its real and imaginary parts for analyzing the dynamics when the motion is 
non-circular. For numerical calculations it is convenient to take all the variables to be dimensionless, so we normalize 
as ip — (j)/m, k — k/m, t — rat, and £ = mx. Writing tp — (ipi + i<p-i)l\/2, we get the equations for the homogeneous 
mode as 



<Pi + 



1 + K + Klog 



P\ 



2M 2 



0. 



where h = H/m, i = 1,2, and the prime denotes the derivative with respect to t, and, for the fluctuations, 



dr 2 



3h T- + "a + V H 
dr a 1 



( 8<Pi 



0, 



(27) 



(28) 



where Vij denotes the second derivative with respect to ipi and tpj, and explicitly written as 



Vh = l + K + K\og 
V 12 = V 2 i = 2K^1 



p\ + pI 

2M 2 



2K 



Pi 



2 1 9 ' 

P1 + P2 



9 ' 

P1 + P2 

In particular, when the homogeneous part is just oscillating along the radial (i.e., real) direction, we have 

~k 2 



Sip" + 3h5<p'i 



1 + /3iK + K log 



2M 2 



5<Pi = 0, 



(29) 



(30) 



where j3i = 3 and 02 = 1. At first glance, we could not see any crucial difference between evolutions of real 
and imaginary parts of fluctuations. As we see in Fig. |l|, however, the imaginary part does not grow at all, while 
considerable growth can be seen in the real part. (Here we have neglected cosmological expansion for simplicity.) 
Notice that both real and imaginary parts of the fluctuations develop when the homogeneous part is rotating in the 
potential.^] (See Fig. ||.) In either case, the instability band coincides with the analytical estimate (26). 

One might thus be tempted to conclude that there is no Q-ball formation, but, as we will see in the next Section, 
mode-mode (self) couplings play a crucial role for developing the imaginary part of the fluctuations, resulting in Q 
and anti-Q balls to form. 



IV. MODE-MODE COUPLINGS 



It is not apparent that the fluctuations in the imaginary direction grows on the top of the pure oscillating homoge- 
neous condensate in the real direction. As mentioned in the previous Section, however, self-couplings of the inflaton 
field are expected to create fluctuations in the imaginary direction as well. This can be best seen in lattice simulations, 
which we have done in one, two, and three spatial dimensions. Since, all the cases are more or less the same in the 
qualitative sense, and the resolution of the three dimensional lattices is not so good, we have grounds for believing 
that the two-dimensional results presented here will also capture the qualitative behavior of the full three-dimensional 
case. In the context of Affleck-Dine baryogenesis, lattice simulations have been done previously |T^,[l^,|l3|,^2],^3) , but 
it is usually assumed that the homogeneous condensate is rotating, not just oscillating in the potential. In addition, 
here the universe is dominated by the inflaton field in contrast to the Affleck-Dine case, in which a matter-dominated 
universe is usually assumed. We also note that Q-ball formation in the pure oscillation of the homogeneous condensate 
has been numerically verified previously in the context of Affleck-Dine mechanism [Q , and discussed in the context 
of quintessence and oscillating inflation pa] . 



3 When the field is rotating in the potential, fluctuations do not grow for positive K. However, the imaginary part of the 
fluctuations actually grows when the homogeneous field is just oscillating. The growth rate is much larger than that for the 
negative K, but, needless to say, Q-ball formation never occurs, since a Q-ball solution does not exist for positive K. 
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FIG. 1. Instability bands when the homogeneous mode is just oscillating along the real direction in the potential (|l|). Solid 
and dashed lines denote the fluctuations in the real and imaginary directions, respectively. Here we take K = —0.02. 
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FIG. 2. Instability band when the homogeneous mode is rotating circularly in the potential ([j]). Solid and dashed lines 
denote the fluctuations in the real and imaginary directions, respectively. Here we take K = —0.02. 



After inflation the inflaton field is almost homogeneous beyond the horizon, and its motion in phase direction has 
completely frozen. However, there exist small variations, which are the seeds for the large scale structures in the 
universe. We can thus pose the initial condition at each lattice site as follows: 



MO) 

M0) 



90 ' 
A 3 , 



Ai, 0i(O) = o + A 2 
02(0) = A 4 , 



(31) 



where tpo = Mp/^Vlix, 0o = — mMp/ \Zl2n, and A's are small numbers which represent small fluctuations. The 
equations to be integrated are 



1 + K + K log 



™ (n + n) 



0. 



(32) 



where i = 1 and 2, which denote real and imaginary parts of the field, respectively, and rh — m/M. Here we will mainly 
discuss results from two-dimensional lattices with the box size N = 2048 and the lattice spacing A£ = 0.004, and 
K = —0.1, and show qualitative three-dimensional results with the box size N = 128, the lattice spacing A£ = 0.01, 
and K = —0.1 as well. 

Fluctuations in the real direction grow very fast, and their amplitudes become of the same order as that of homo- 
geneous mode at r ~ 1000. It is depicted in Fig. [| where we show the time evolution of both real and imaginary 
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FIG. 4. Power spectra of fluctuations in the real and imaginary directions just before the real part fully develops. 



parts of the fluctuations. Fluctuations in the imaginary direction, however, have not yet grown at that time, and it 
starts growing later. 

In Fig. [|, we plot the spectra of real and imaginary parts of the fluctuations at t = 826, just before the former fully 
develops to the size of the background field value. In fact, the instability band lies in the right range k < a(r)|i(r| 1 / 2 TO, 
and the most amplified mode coincides with the anticipated Q-ball size. After the real part of the fluctuations fully 
developed, other modes outside the instability band also grow due to mode-mode couplings, as can be seen in Fig. |g, 
which shows the spectra of the fluctuations at r = 1098, just after the real part has fully developed. Although the 
spectrum becomes wider and smoother, the position of the most amplified mode, which corresponds the Q-ball size, 
stays still. 

One surprising fact is that lumps are actually formed just after the real part of the fluctuations has developed fully 
even though the imaginary part has not yet started to grow. But, these lumps have not yet contained charges inside. 
This can be gathered in Figs. ^ and fj], where we respectively show the charge density and the amplitude of the field 
at the time just after the amplitude of the fluctuations in the real direction catches up with that of the homogeneous 
mode. Although the field has not settled down completely, we can already see many lumpy objects, which have no 
charge. Notice that they have in fact the right size of the prospective Q balls: ~ \K\~ 1 / 2 m~ 1 . 

As time goes on, charges are accumulated into lumps to become Q balls due to mode-mode couplings between 
fluctuations in the real and imaginary directions. We show the evolution of (<5cj 2 ) in Fig. ||. It is evaluated from 
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FIG. 6. Charge density of the field. Here we show only a small portion of the entire lattice of 8.192 x 8.192. 
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FIG. 7. Amplitude of the field with the same range as in Fig. 
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FIG. 8. Evolution of the fluctuation of angular velocity. 
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FIG. 9. Power spectra of fluctuations in the real and imaginary directions just after the real part fully developed. 



(fe 2 ) 



-<^i> 



^2 



(M 2 



(33) 



We see that (8uj 2 ) starts to grow right after the field fluctuations in the real part have fully developed. We may regard 
it as a local evolution of the charge distribution, or, in other words, the charge accumulation into lumps, which leads 
to the Q-ball formation. Once it is saturated (at r ~ 3500), lumpy objects will not disintegrate but remain as Q balls 
because of the charge conservation. We can also see in Fig. ^, where the spectra of fluctuations in real and imaginary 
parts are shown. At around this time, the amplitude of the imaginary part has almost matched with the real partner, 
which also implies the actual Q-ball formation. 

At this time, the field values settle down, and Q balls and anti-Q balls are completely isolated from the background. 
Figures [lO] and |ll| show the charge density and the amplitude of the field after the imaginary part of the field fluctuation 
has fully developed. Now we can see "mature" Q balls, and charges of the Q balls conserve from that time on. 

Of course, a Q-ball formation is seen on the three-dimensional lattices. We only show how it appeared in Fig. [T^. 
Here, we plot the configuration at r = 3860, just after Q balls are completely formed. The typical size and charge of 
the Q balls are R p hy S / m ~ 0-05 x a(3860) ~ 10, which almost coincides with the analytically estimated values, and 
Q ~ (a few) x 10 6 

Therefore, it appears that in running mass inflation Q balls are produced not only by the negative pressure which 
makes the homogeneous condensate to fluctuate at the same scale as the Q-ball size, but the mode-mode coupling 
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FIG. 10. Charge density distribution in a small sub-lattice. 



t=3546 





FIG. 12. Q balls in three-dimensional lattices with N = 128 and A£ = 0.01 at r = 3860. 
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plays a crucial role in creating locally distributed charge density. In this sense we can state that in a potential such 
as Eq.(|l|) a kind of preheating leads to Q-ball formation. 

V. REHEATING THROUGH EVAPORATION 

So far, we have seen it is very natural for the inflaton condensate to deform into Q balls in a certain inflaton 
potential such as Eq.(Q). As a consequence, it may provide the alternative way of reheating after inflation, namely 
reheating through evaporation, which leads to slower reheating process than usually considered, as we claimed in 
Ref. [Q. Q balls in the potential (]]]) are unstable for decaying into particles lighter than ~ m Jl5fl . If the decay 
products are pairs of fermions, the decay process occurs only on the surface of Q balls because of the Pauli blocking 
inside the Q ball @. 

A Q ball of size R ~ \K\~ 1 l 2 m~ 1 forms when the fluctuations grow nonlinear. Since the growth rate of fluctuation 
is ~ IA'Ito, the Hubble parameter at the formation time will be estimated as 

H f ~ j\K\m, (34) 

where 7 ~ 2 x 10~ 2 is some numerical factor derived from numerical simulations. As we assume that \K\ <C 1, we 
can approximate the decrease in the amplitude of the oscillations by <pf ~ 4>i{H f / Hi) as in the matter dominated 
era, where fa ~ Mp/^/V2n, denotes the amplitude at the end of inflation in chaotic model, and Hi ~ to when the 
oscillations begin. Therefore, the total charge of a Q ball is given by 

Q „ ^i?3 n? „ lp?f\K\*I&mM*, (35) 

where n q = ftuxpQ, (fto ~ C4>f, and (3 -C 1 and C, 1 are numerical factors. 

As we shall see below, the evaporating reheat takes over the usual perturbative decay in relatively large coupling 
where h < 1. In this regime, the evaporation rate is saturated by 

dQ uj 3 A 

dF-192^' (36) 

where A ~ 4nR 2 is the surface area of the Q ball, and w ~ m. Thus, we obtain the decay rate of the Q ball through 
evaporation as 

_l_dQ 3 / m \ 

Q Q dt ~ 16^/3CVI^I 3/2 \Mp) m ' { ' 

Note that the decay rate is determined by the ratio m/Mp ~ 10 -6 , which is fixed by the anisotropics seen in the 
cosmic microwave background radiation. Therefore, even though we are in relatively large coupling limit, the decay 
rate mimics that of a Planck suppressed interaction of the inflatonic Q ball to the matter fields. 

Now let us compare the rate (|37|) with the perturbative decay rate T pert = h 2 m/8Tr. Taking the ratio, we have 



Tq 3 / TO x 2 



r pert 2PC 2 l 2 \K\ 3 / 2 h 2 \M P 



(38) 



When this ratio is less than unity, the evaporation process from the Q-ball surface suppresses the reheating. This is 
achieved if h > 10 -3 for reasonable parameters, such as \K\ ~ 0.1, f3 ~ 10 -2 , 7 ~ 2 x 10~ 2 , and £ ~ 2. These values 
of the parameters are derived from numerical simulations on the three-dimensional lattices: Q ~ 3 x 10 6 , for example. 
Since the reheating temperature can be estimated as 



T RH ~ Q.ly/TMp oc T 1 / 2 , (39) 

it will be suppressed by a factor ~ 2 x 10 -3 for h ~ 1. 

In general K and h are not independent quantities but are related to each other by \K\ ~ Ch 2 /16ir 2 , where C is 
effective number of bosonic and fermionic loops. If the inflaton sector does not belong to a hidden sector, it is very 
natural that the inflaton coupling to other matter fields is relatively large, i.e. h > (m/Mp). We can thus have easily 
a situation where \K\ ~ 0.1, if, e.g., C ~ 10 and h ~ 1. 

To end this section, we will comment on (fermionic) preheating pq , ^7| . In the preheating stage, it is very easy to 
transport inflaton energy into a pair of massless fermions in an efficient way with the maximum physical momentum 
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k ma x cx {Wm) 1 ' 2 , (40) 

where 4> is the amplitude of the oscillating inflaton. Even though fcrmionic preheating is efficient, the whole inflaton 
energy is not transferred in this process. The energy density stored in the fermions remains small compared to the 
inflaton energy density, i.e. the comoving number density and the energy density of the massless fermions follow 

k 3 max {h^mf 2 ki ax h 2 <j> 2 m 2 

nf ~-&^~—^z—' Pf ~-8^~~8^- (41) 

The ratio of fermion energy density compared to the inflaton energy density after the end of fermion preheating 
behaves as 

(42) 

P 87T Z 

Hence most of the energy density of the inflaton does not transfer to that of the fermions. This result is quite robust 
and remains valid even if the fermions have a non-zero bare mass p8| . 

Even if fermions are produced, they cannot scatter inflaton quanta off its condensate p6[ , unlike in the case of 
bosonic preheating with an interaction g 2 (f> 2 x 2 - The zero mode of the inflaton thus remains intact as a coherent 
condensate which will fragment into Q balls later. Notice that, as mentioned earlier, a Q-ball formation is reminiscent 
of bosonic preheating due to the presence of attractive self coupling of the inflaton, which stems in the logarithmic 
term in the potential. Therefore, Q-ball formation appears to be a robust feature of any inflationary model with a 
potential allowing for a Q-ball solution. Moreover, for a flatter potential, the size of the Q ball is correlated with the 
charge, as in the context of the Affleck-Dine condensate with the gauge-mediated SUSY breaking. This means that 
the suppression on the reheating will be much more efficient in those cases. 



VI. CONCLUSION AND DISCUSSION 



Inflaton condensate fragmentation and the formation of lumps may be a rather generic, albeit in many cases a 
transient feature of preheating. However, for the running mass inflaton with K < 0, and possibly for other types of 
inflation admitting at least temporarily a mass term growing slower than |$| 2 , the final outcome is markedly different 
from the usual preheating. The effective equation of state, averaged over field oscillation in the real direction, has 
negative pressure, which is the cause of the condensate break-up. If the inflaton happens to be a complex field so 
that it admits soliton solutions carrying a global charge, one would typically expect an inflatonic Q and anti-Q ball 
pair to be generated when within a Q-ball forming region the field fluctuations in both the real and the imaginary 
directions exceed the coherent background field value. 

To verify the generation of inflatonic solitons it is not sufficient to follow only the evolution of the field amplitudes. 
In the real direction the fluctuations grow very fast to the size of the prospective soliton lumps, but this by itself is 
not sufficient to ensure Q-ball formation. Indeed, field fluctuations in the imaginary direction do not grow at all in 
the linear approximations. To see what really happens, one has to perform simulations in a configuration space on 
the lattices, as we have done in the present paper. 

In our simulations we can see the formation of lumpy objects but already at relatively early stages when the mean 
fluctuations in the imaginary direction have not yet develope. These objects have not become Q balls yet, but they 
eventually accumulate charges to become Q and anti-Q ball pairs. We believe that the reason is that there will be 
large local fluctuations which, once pushed over the limit after which bound solitonic states can exist, cannot fluctuate 
back. The attractive inflaton self-interaction generates locally distributed U(l) charges through mode-mode couplings 
which then evolve into semi-stable inflatonic Q balls. The formation of such solitonic states is non-trivial and stands in 
contrast to the lumps of the K > case; the latter is likely to be an evanescent feature that would merely disintegrate 
and decay away very fast. 

Once produced, these Q balls will decay by evaporation and hence give rise to a "slow" reheating because of the 
suppression which comes from the surface-volume ratio. In particular, this effect becomes more efficient for rather 
larger coupling h ~ 1, which might help the disastrous situations such as the gravitino problem in the supersymmetric 
scenarios. In fact, the surface evaporation through Q ball is perhaps the most natural way to obtain low reheat 
temperatures inspite of having order one inflaton coupling to the fermions, which solves not only the gravitino over- 
abundance problem but also moduli, and dilaton, especially when the inflationary scale is sufficiently large. 

Although we have looked upon the running-mass potential, existence of the Q-ball solution is a very generic feature 
in a certain class of models. It would be interesting to study whether such inflatonic Q balls could be seen to arise also 
in other inflationary models such as hybrid inflation, new inflation, etc., and study in these directions are in progress. 
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